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2 ■ ABSTRACT 

^ I We report the identification of SH A?; = 1 ro-vibrational lines in the pubhshed high- 

ly ' resolution infrared spectrum of the S-type star, R And. This is the first astronomical 

■ detection of this molecule. The lines show inverse P-Cygni profiles, indicating infall 

■ motion of the molecular layer due to stellar pulsation. A simple spherical shell model 
with a constant infall velocity is adopted to determine the condition of the layer. It 
is found that a single excitation temperature of 2200 K reproduces the observed line 
intensities satisfactory. SH is located in a layer from 1.0 to ~ 1.1 stellar radii, which 
is moving inward with a velocity of 9 kms~^. These results are consistent with the 
previous measurements of CO Au = 3 transitions. The estimated molecular abundance 
SH/H is 1 X 10"'', consistent with a thermal equilibrium calculation. 

Subject headings: infrared: stars — line: identification — stars: AGB and post-AGB — 
stars:atmospheres — stars: individual(R And) 
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1. Introduction 

Chemistry of sulfur-containing species in space is especially interesting due to its chemical 
activity and relatively high abundance (Duley, Millar, & Williams 1980). It is known that abun- 
dances of sulfur-bearing molecules are strongly influenced by the presence of shocks (Hartquist, 

Oppenhcimcr, & Dalgarno 1980). So far, 14 sulfur-bearing species have been identified in space; 
CS, SO, NS, SiS, S0+, SH2, OCS, CCS, C3S, H2CS, HNCS, SO2, HCS+, and CH3SH. However, 
the simplest sulfur compound, SH, has not been detected in spite of radioastronomical searches 
using the A-type doubling transitions (Meeks, Gordon, & Litvak 1969; Heiles & Turner 1971). 

While sulfur chemistry in the circumstellar envelope of red giant stars has been studied ex- 
tensively (e.g., Omont et al. 1993), that in the atmosphere has been regarded as in the thermal 
equilibrium state (Tsuji 1964, 1973). In the atmosphere of oxygen-rich giants, SH is the molecule 
first formed from the sulfur atom, at the temperature of about 2000 K. The SH abundance decreases 
when temperature is below ~ 1500 K, and H2S and SiS become the dominant species. On the other 
hand, Yamamura et al. (1999a) report the detection of the 1^3 infrared band of SO2 in the spectra of 
oxygen-rich Mira variables obtained by the Short- Wavelength Spectrometer (SWS: de Graauw et 
al. 1996) on board the Infrared Space Observatory (ISO: Kessler et al. 1996). Their model analysis 
indicates that the molecules are located in the extended atmosphere, at about five stellar radii. The 
abundance of SO2 is 5-10 orders of magnitude larger than the values in thermal equilibrium (Tsuji 
1973; Woitke et al. 1999). The result strongly STiggcsts the presence of non-equilibrium processes 
in the extended atmosphere, probably related to shocks due to stellar pulsation (Beck et al. 1992; 
Duari, Cherchneff, k Willacy 1999). 

In this paper, we report the detection of SH ro-vibrational transition lines in the published 
high-resolution infrared spectrum of the S-type star, R And. We discuss the physical and chemical 
conditions of the SH layer in the star. 

2. Identification of SH lines 

The data used in this study were obtained with the Fourier-transform spectrometer at the 
Kitt Peak National Observatory 4 m telescope, and were published in Ridgway et al. (1984). R 
And is a Mira variable with a period of 409 days and a visual amplitude of 9.1 mag (Kholopov 
et al. 1988). The spectral type ranges in S3,5e-S8,8e. The observation of R And was performed 
at the optical variable phase of (p = 0.82. Unfortunately, the reduced data were lost after several 
updates of archive formats at the observatory. Therefore, we re-digitized the data from the printed 
spectrum. The results may be as accurate as the original within ~ 0.02 cm~^ in wavenumber and 
a few % in intensity. The spectrum was shifted by 5.2 kms~^ to longer wavelengths to adjust the 
absorption minimum of HCl lines at their rest frequencies (Ridgway et al. 1984). No correction for 
terrestrial motion was applied. In their original paper, Ridgway et al. (1984) give identifications 
of OH, NH, CH, SiO, CS, HCl, and atomic lines. The spectrum of R And, especially in the 
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wavelengths shorter than the SiO first-overtone bandheads, is dominated by HCl fines. OH and 

NH fines are detected, but much less prominently than in the oxygen-rich stars like a Ori and a 
Tau. There arc many strong lines left unidentified. Figure 1 sfiows tfie spectrum of R And between 
2700 and 2750 cm~^. The positions of SH and HCl ro- vibrational transitions are indicated. The 
frequencies are calculated from the molecular constants given by Ram et al. (1995) for SH, and are 
taken from HITRAN database 1996 edition (Rothman et al. 1992) for HCl. It is obvious that many 
unidentified strong lines are attributed to SH transitions. We assigned 39, 41, and 11 transitions 
in the v = 1 — 0, 2 — 1, and 3 — 2 bands, respectively, between 2500 and 2778 cm^^. The lower 
state energies of these transitions are up to 7000 cm"-*^, indicating that the molecules are highly 
excited. The isotope ^^SH may also be detected. The frequencies of ^^SH lines are estimated from 
the molecular constants of SH using the relation of reduced mass (Herzberg 1950). 



3. Modeling 

We analyze the observed SH line profiles with a simple model. The SH (and also HCl) lines 
in R And show inverse P-Cygni profiles, indicating that the molecular layer is moving inward to 

the star due to stellar pulsation. We, therefore, apply a spherical shell model with a constant infall 
velocity. The molecules in front of the star cause absorption, while those extended in the blank 
sky contribute as emission. Considering that the shell may be in the atmosphere, wc adopt an 
exponential density law, n(r) oc exp(— r/ro), where ro is the scale height. The star is assumed 
to be a 3000 K blackbody. For simplification, a constant excitation temperature is adopted, and 
the energy population of the molecule is calculated by assuming local thermodynamic equilibrium 
(LTE). This may be justified if the molecular shell is thin and in the high density region near the 
photosphere. The line intensity of SH is calculated based on Benidar et al. (1991), which takes 
account of the Herman- Wallis effect. The spectrum is normalized by the stellar continuum. No 
smoothing is applied. 

The excitation temperature can be determined from the relative intensities of the lines at 
different energy levels; higher temperatures excite the molecules to higher energy levels and increase 
the line strength from these levels. We find that 2200 K is most reasonable for the present case. 
A turbulent velocity of 6 kms~^ and infall velocity of 9 kms~^ reproduce the line profiles. The 
model spectrum is shifted by —14 kms~^. Considering the terrestrial velocity of —2.7 kms~^ at 
transit of the observation, and the shift applied by Ridgway et al. (1984), +5.2 kms~^, a radial 
velocity of —16.5 kms^^ is obtained. The inner radius and the scale height of the SH shell is 1.0 
and 0.08 times the stellar radius, respectively. The SH column density in the shell is 4.0 x 10^° 
cm~^. The uncertainties of the parameters, which changes the model spectrum by about 10 %, 
are ±100 K for excitation temperature, a factor of two for column density, ±0.01 R^ for scale 
hight, and ±1 kms^^ for the velocities. The radial velocity has an error of ±2 kms~^ due to the 
re-digitization process. The model spectrum is compared with the observation in Figure 2. The 
wavelength regions are selected so that the contamination of other spectral lines is minimal, and 
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that a wide range of energy levels is covered. Strong lines in Figure 2 are listed in Table 1. The 
fit is satisfactory in most of the lines. Fitting of "'^SH lines results in an isotopic abundance ratio 
of 5~10 %. We note that the same parameters also give reasonable fits for the HCl lines with a 
column density of 3.0 x 10^^ cm~^. 

4. Discussion 

Why are SH lines so prominent in R And, an S-type star? S-type stars are characterized by 
their chemical anomaly of similar carbon and oxygen abundances. This leaves few C and O atoms 
for further chemical processes, after formation of CO molecules at 4-5000 K. SH should be quite 
abundant in the atmosphere of oxygen-rich stars, but a minor product in a carbon-rich environment 
(Tsuji 1973). Since we do see SiO lines, but no CS line, in the spectrum of R And, the star is 
slightly oxygen-rich. Nevertheless, it is not expected that SH is much more abundant in S-type 
stars than 0-rich stars. This also holds for HCl, which is the most abundant Cl-bearing molecule in 
an oxygen-rich environment, although the absolute abundance is by one order of magnitude lower 
than SH. Ridgway et al. (1984) suggested that the atmosphere of S-type stars is more transparent 
than 0-rich stars, so that the lines of minor species could be stronger due to larger path length. 
We especially emphasize the contribution of the H2O molecules. In the same paper, Ridgway et 
al. reported that an enormous number of lines heavily blanket the 2400 2800 cm^"*^ region of the o 
Cet spectrum. They suspected that these are highly excited water lines. This is supported by the 
ISO/SWS observation of o Cet (Yamamura et al. 1999b). Despite the relatively poor resolution of 
the SWS spectrum, they demonstrate that the star shows highly excited H2O lines in the 3.5-4.0 
^m region, arising from hot (2000 K) and optically thick molecular layer. Probably, in 0-rich stars, 
SH (and HCl) lines are weaker because of the shorter path length, and also heavily contaminated 
by the hot water lines. On the other hand, H2O molecules are not favored in the atmosphere of 
S-type stars. The thermal equilibrium abundance of H2O in S-type stars is 3-4 orders lower than 
in O-rich stars (Tsuji 1964). 

The derived infall velocity and the excitation temperature of the SH layer are consistent with 
the measurement of CO At> = 3 lines by Hinkle, Scharlach, h Hall (1984) at ^ = 0.80, indicating 
that the SH molecules are in the same region as hot CO near the photosphere. Comparison with 
the CO column density, 2.6 x 10^^ cm~^, given by Hinkle, Scharlach, & Hall (1984), leads to 
SH/H = 1 X 10"''. In thermal equilibrium at the temperature of 2200 K and gas pressure of 
log(Pg) = 1.0, the SH abundance is 2 x 10~^ (Tsuji 1973, 1999). This value is insensitive to the 
C/0 ratio as long as C/0 < 1. The present estimate of SH abundance in R And, although it is 
rather crude, is consistent with this calculated value. The sulfur chemistry in SH layer still follows 
thermal equilibrium, because of its high density and high temperature. 

We see no clear evidence of distinct velocity and/or different temperature component of SH 
lines in the spectrum. This implies that the SH molecules are distributed only in a thin layer in 
the atmosphere. In thermal equilibrium, SH is most abundant around 1800 K, and then is rather 
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quickly transformed to H2S or SiS below ~ 1500 K (Tsuji 1973). Otherwise, non-equilibrium 
chemical reactions may lead to completely different compositions in the extended atmosphere, e.g. 
the enhancement of SO2 in oxygen-rich stars (Yamamura et al. 1999a; Beck ct al. 1992). We could 
not find any clear indications of H2S or SO2 in the spectrum of R And. The upper limit of the H2S 
column density is at least a factor of 10 larger than the observed amount of SH. A dioxide molecule 
SO2 may not be abundant in the atmospheres of S-type stars. Other possible candidates, SO or 
SiS, have no transition in the present wavelength coverage. 

The authors arc grateful to Prof. T. Tsuji for his suggestion on thermal equilibrium chemistry 
and calculations of SH abundance. I.Y. acknowledges financial support from a NWO PIONIER 
grant. 
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Fig. 1. — A part of the spectrum of R And obtained by Ridgway et al. (1984) is presented. The 
positions of Av = 1 transitions of SH and ^^SH (upper), and HCl (lower) are indicated by ticks. 
The notation of the transitions is the electronic angular momentum, the upper vibrational level, 
and the rotational line index for SH, and the rotational line index, the isotope index (5: H^^Cl, 7: 
H^^Cl), and the upper vibrational level for HCl, respectively. 

Fig. 2. — The synthesized model spectrum (upper) is compared with the observation in six wave- 
length regions. The regions are chosen so that contamination of other spectral lines is minimum, 
and that a wide range of energy levels is covered. The model assumes a spherical shell falling toward 
the stellar surface with a velocity of 9 kms~^. The molecules are distributed from 1.0 R* with a 
scale height of 0.08 R*. The model assumes constant excitation temperature of 2200 K, column 
density of 4.0 x 10^'^ cm~^, and Dopplcr broadening velocity of 6 kms^^. The model spectrum is 
shifted by —14 kms~^. The lines indicated by the ticks are listed in Table 1. 
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Table 1. SH lines indicated in Figure 2. 



No.^ Frequency Transition'^ Ei^ 

cm"~^ cm~^ 



1 


2532, 


.1257 


3/2 


1 


— 


P 3.5ff 


110, 


,7175 


1 


2532, 


.I486 


3/2 


1 


— 


P 3.5ee 


110. 


.6810 


2 


2534, 


.7858 


1/2 


3 


— 2 


R 7 5ee 


6038. 


,9539 


2 


2534, 


.8657 


1/2 


3 


_ 2 


R 7.5fF 


6040. 


,7705 


3 


2546, 


.2628 


1/2 


2 


- 1 


R 1 5ee 


2994. 


.8376 


4 


2546, 


.4908 


1/2 


2 


— 1 


R 1.5fF 


2995. 


,3755 


5 


2546, 


.9511 


1/2 


3 


_ 2 


R 8 5ee 


6192. 


,1858 


5 


2547, 


.0067 


1/2 

X/ z, 


3 


_ 2 


R 8.5ff 


6194. 


,1577 


6 


2549, 


.2843 


1/2 


1 


— 


P 2.5ff 


446, 


.0913 


7 


2549, 


.5750 


1/2 


1 


— 


P 2.5ee 


445, 


.2626 


8 


2641, 


.3159 


3/2 


3 


_ 2 


E24 5ee 


10354. 


,2802 


8 


2641, 


.3458 


3/2 


3 


— 2 


R24.5fF 


10358. 


,7197 


9 


2642, 


.8296 


3/2 


1 


- 


E 1 See 


0, 


.0000 


9 


2642, 


.8395 


3/2 


1 


— 


R 1.5ff 


0, 


.0037 


10 


2644, 


.6568 


1/2 


1 


— 


R 1 5ee 


396, 


.9951 


10 


2644, 


.7048 


3/2 


2 


— 1 


R 8 5ee 


3286. 


.0606 


10 


2644, 


.7949 


3/2 


2 


— 1 


R 8.5ff 


3286, 


4309 


10 


2644, 


.8958 


1/2 


1 


-0 


R 1.5ff 


397, 


.5540 


11 


2659, 


.3953 


3/2 


1 


-0 


R 2.5ee 


46, 


.1293 


11 


2659, 


.4143 


3/2 


1 


-0 


R 2.5ff 


46, 


.1440 


12 


2661, 


.5489 


1/2 


2 


- 1 


R 9.5ee 


3889. 


.9056 


12 


2661, 


.5866 


1/2 


2 


- 1 


R 9.5fF 


3892. 


.0848 


13 


2661, 


.9331 


1/2 


1 


- 


R 2.5ee 


445, 


.2626 


13 


2662, 


.1518 


1/2 


1 


- 


R 2.5ff 


446, 


.0913 


14 


2709, 


.9111 


1/2 


1 


- 


R 5.5ee 


705, 


.7060 


14 


2710, 


.0575 


1/2 


1 


-0 


R 5.5ff 


707, 


.2638 


15 


2710, 


.9529 


1/2 


2 


- 1 


R14.5ff 


5042. 


.3624 


15 


2711, 


.0189 


1/2 


2 


- 1 


R14.5ee 


5039. 


.8647 


16 


2770, 


.9408 


3/2 


1 


- 


R10.5ee 


1074. 


,9261 


16 


2771, 


.0636 


3/2 


1 


-0 


R10.5ff 


1075. 


.6188 



'^Numbers in Figure 2. 

^Electronic angular momentum index. 

'^Vibrational and rotational transition indices. 
"^Lower state energy relative to v = 0,J = 1.5,0 = 



